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Some analytical solutions for spherical-symmetrical equilibrium configu- 
rations (planets/stars) in Newtonian Theory of Gravitation (NTG) with de- 
viations from 1/r law are discussed. Stability of star against the first-order 
phase transition is particularly influenced by deviations of degree in power-law 
potential . 
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I. INTRODUCTION 

In the classical NTG, with the point-mass potential 

U(r) = Gm/r, (1) 



^ | the equation of the spherical-symmetric equilibrium states reads: 



1 dP{r) Gm{r) 

(2) 



Q-f p(r) dr r 2 
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Here, as usual, G stands for constant of NTG, P(r), p(r) are local pressure and density of 
stellar matter, at distance r from the center of star, and m(r) is mass inside the sphere of 
radius r ("r-sphere"). Provided EOS P(p) is known, the Eq. (2) is ODE and is degenerate 
in the sense, that any deviation from standard law (1) spoils the fundamental characteristic 
of spherical-symmetric stars - dependence of local gravitational force, at radius r, only on 
the mass inside "r-sphere". That is any deviation of 1/r law leads, in principle, to the 
enormous complexity of problem of the spherically-symmetric equilibrium states. Instead of 
ODE one has to deal with integro-differential equation (IDE), see next section. 
In [1] authors discuss some recent schemes of compactivation leading to the generalized form 
of the point-mass potential, 

tf(r) = — (l + ae-'A), (3) 
r 

where a and A are constants. In the celestial mechanics, the various "non- Newtonian" or 
"non-gravitational" forces have been considered quite a long ago [2]. 

In this note I consider only the problem of stellar structure in NNG (Non-Newtonian Grav- 
itation). As nobody knows the exact form of "non- Newtonian" law, it is reasonable to 
consider first the various simple cases of potential U(r) = Gm * f(r) . Here I consider 
"power-law case", 
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f(r) = 1 / r 1+ «, (4) 

where a is constant. Some manifestations of this law in stellar structure are rather dramatic 
that may help to judge in favor or against it. Despite this "applicational" aspect it's cer- 
tainly worth by itself considering the problem of stellar structure with non-l/r gravitational 
potential. 

At first, I present the hydrostatic spherical-symmetric equilibrium equation in general case. 



II. EQUATIONS OF STELLAR STRUCTURE IN NNG 

If f(x) is "generalized Newtonian" law, then the potential of thin spherical shell is as 
follows: 



u s heii(i",a) da = 2 % G p(a) a 2 da / dx. (5) 

J\r-a\ r a 



Here, a and da are radius and "thickness" of infinitesimally thin shell, p(a) is matter density 
at radius a, and r is distance from the center of spherical shell. 

Note, that only at f(x) = 1/x 1 (and, of course at f(x) = l/x°, when point-mass potential 
itself is constant!), the potential inside the spherical shell is constant (and gravitational force 
is zero). It's worth mentioning that both laws (3) and (4) give the first-order correction 
leading to necessity of IDE, while the zeroth-order term coincides with 1/r law. 
The potential U sp here(p, r, R) , at distance r, of a sphere of radius R is : 

U sp here(p, r, R) = I u sheU (r, a)da. (6) 

JO 

The Eq. (6) is valid both inside and outside the sphere. 

The equation of the hydrostatic equilibrium of spherical-symmetric star reads: 



1 dP(r) _ dU sphere (p(r), r, R) 



p{r) dr dr 



(7) 



The right side of this equation is a double integral, and in general is of much more compli- 
cated form than in the standard f(r) = 1/r case, Eq. (2). Provided EOS P(p) is known, 
Eq. (7) is IDE for density function p(r). 

As at any arbitrary EOS P(p) and non-l/r law it's apparently impossible to solve the prob- 
lem (5-7) , it's a good idea to consider some simple cases first. In the next section, I present 
the spherically-symmetric model of a homogeneous imcompressible liquid, that corresponds 
to a polytropic index n = in the polytropic EOS of the form P(p) = K p 1+n . 

III. HOMOGENEOUS SPHERE 

At the point-mass potential of form (4), the gravitational potential of homogeneous 
sphere of radius R and density p = const is as follows: 

aN = N-a; t = r/R; 
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(8) 



Here upper and lower signs correspond to t < 1 and t > 1 respectively. Values of potential 
at the center and the surface of the homogeneous sphere are: 

U(t = 0) = AnGpR a2 /a2] U(t = 1) = 27rGp(2R) a2 /(o2 a3). (9) 

The total gravitational potential energy of homogeneous sphere is: 

W = -7r 2 G'p 2 (2 J R) a5 /(a2 a3 a5). (10) 

Solution of equation (7) in this case leads to the following relation for central pressure: 

(11) 



P c = p* [U(t = 0) - U(t = 1)] = 2nGp 2 R a2 (2 - 2 a2 /a3)/a2. 
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FIG. 1. Dependence of the potential (left panel) and the gravitational force (right panel) inside 
the homogeneous sphere for \/r 1+a point-mass potential. Vertical scales are somewhat arbitrary 
as for different values of a dimensionality of constants are different, see Eqs. (8-11). 

In Fig. 1, some relevant curves are presented. Of course, in "real" cases it appears that 
|a| << 1, so the Fig. 1 is presented here only for pure illustrative purposes. We note that 
nothing extraordinary interesting happens with the gravitational potential of homogeneous 
sphere by broken 1/r-law even for such large values of \a\. Suprisingly enough, even a such 
particular characteristic of NTG as maximal force at surface of sphere survives in NNG. 
It may be noted, that in spite of "enormous", in principle, complexity of the spherical- 
symmetric equilibrium equations in NNG, in reality it is difficult to point out any macro- 
effect of "micro" -values of a. 

By the way, we notice a rather common opinion that in astronomical ("cosmic") scales, even 
small deviations from 1/r law may lead to some macro effects. For example, external parts 
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(outside the "i?-sphere") of homogeneous dark-matter or any other background matter are 
usually assumed as non-influencing the equilibrium and stability of spherical structure with 
given radius R [3]. 

It's seems worth reconsidering this problem with account of deviations from 1/r law. Two 
potentially important cases may be the models with adiabatic index of EOS close to 4/3, 
suggested sometimes as model of supermassive stars (quasars), and the isotermic configura- 
tions suggested as model of inter-(or proto-) stellar clouds or other cosmic structures [9]. 
Here I consider one rather interesting effect of NNG. In NTG (as well as in GR) the problem 
of equilibrium and stability of a star (planet) with first-order phase transition (PT1) exists 
apparently since pioneer works by Ramsey [4]. 

In "classical" NTG, the elegant and quite suprising result is that if PT1 occurs at the star's 
center with q > 3/2, q = p 2 /pi being the ratio of new-to-old phase densities, the star un- 
dergoes a stability loss irrespective of EOS! 

This "exact" value, 3/2, has apparently its origin from another two exact values: 1, power 
of r in point-mass potential law, and 3, dimension of space. It's the place here to illustrate 
the first part of this assertion, for the simple model. 

The simplest model of star (planet) with PT1 is a two-zone model with constant densities 
Pi and p 2 in envelope and core respectively. 

IV. TWO-ZONE MODEL WITH PT1 

In this model density distribution is a piece-wise function: p(0 < r < r n ) = p 2 = const, 
and p(R > r > r n ) — p 1 — const, with q = p 2 j p\ > 1; here r n and R are radii of core and 
of total star. At the core-envelope boundary, the pressure is constant, P(r n ) = P = const. 
The potential in envelope may be considered as a sum of internal potential of homogeneous 
sphere of density p\ and radius R and of external potential of homogeneous sphere of density 
p 2 — Pi and radius r n : 

Uenv{p\,T, R) = U int , S phere(Pl,r, R) + U out , sp here{p2 ~ Pi, T, r n ), (12) 

where two functions in the right site are given in Eq. (8). Hydrostatic equilibrium equation 
(7) has the solution, in the envelope, r > r n : 

P{r)= Pl [U env ( Pl ,r,R) - U env ( Pl ,R,R)\; (13) 

[Uenv \Pl j T n , 

R)- U env ( Pl ,R,R)}. (14) 

The last expression is "initial condition" for the hydrostatic equilibrium equation (7) and 
determines the total radius R as function of radius of core r n . Total mass of model is: 

M = 4 vr G [ Pl R 3 + (p 2 - pM/3. (15) 

On the other hand, Pq is central pressure for "initial" homogeneous configuration with 
"initial values" of total mass M , total radius R , which are defined as in Eq. (11) with 

R = Rn'. 
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P = 2nGp(R^(2 - 2 a2 /a3)/a2- M = 4 n G Pl B% /3. (16) 

These relations allow to exclude P from Eqs. (13-15) and to express R and r n in units of 
R , and also M in units of M . 

According to the static criterion [5], the dependence of a total mass of equilibrium states 
on central pressure determines the stability of models. That is the branch of M eq (P c ) with 
positive derivative dM eq /dP c presents the stable equilibrium states while the branch with 
negative dM eq /dP c presents unstable equilibrium states. 

In our case we can use r n (or even r n /R) as independent variables as they both are monotonic 
functions of P c . 



A. Classic NTG Case 

For convenience, I present here, very briefly, some relevant formulae (and results) in 
NTG. Although they are known since quite a long ago [3] they have been rederived in 
literature sometimes, e.g. [6]. 

Instead of (13-16) we have (all variables are expressed in relevant "initial" values): 

M — ((q— l)x 3 ) + l)i? 3 , x = r n /R; (17) 

i? 3 - [1 - (2q - 3)r 2 n ]R - 2(q - l)r 3 = 0; (18) 

R=[l + (2q- 3)x 2 - 2(q - l)x 3 }- 1/2 ; (19) 

These equations allow a full analytical treatment of the problem. In particular: 

dM(x)/dx may be negative only at q > 3/2; 

dM(x)/dx = at x determined by the equation 

(q - 1) V + A(q - l)x - (2q - 3) = 0; 

the last equation has solutions x < \/2 — 1; 

at mass minima, M = R 2 . 

B. The l/ r 1+a Case 

As relevant formulae are rather cumbersome in general case I'll not rewrite them from 
Eqs. (8-16), and only consider a situation with small cores. For small values of a relative 
radius of core, x = r n /R « 1, we have expansion up to x 3 : 



R= L. <3-«)(2-«) 



6 - 2a - 2 2 ~ a 



2 2 -°(<7-l) x 2- a %-V X ~ ! f2 "" ) 



(2-a)(3-a) 



(20) 



This equation is valid for arbitrary a and for small x. At the limit a — > this equation coin- 
cides with exact expression (19) which is valid for any x, in NTG case, providing q cri = 3/2. 
However the case of non-zero a's differs from the classic case qualitatively. That is we don't 
have, at small values of a, some small correction to critical value of PT1 in NTG q cr = 3/2. 
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Instead, at a > , the smallest power of x in Eq. (20) is a2 < 2, the corresponding 
coefficient is negative for any q > 1, therefore we have dR/dx < 0. As evident from (15), in 
M "additional" term oc x 3 , that is AM oc AR oc (x 2 ~ a or x 2 ), that provides the continuity 
of derivative dM/dR along the curve of equilibrium states [7]. 

In other words, at a > (more steeper dependence of potential as compared with 1/r law of 
NTG), we have q cr = 1, that is for small cores there is a situation of the absolute instability 
of star against PT1; while at a < 0, there is reverse situation of the absolute stability of 
star against PT1. 

For the larger (but still small enough) dimensions of core there takes place a restoration of 
the normal behavior of M(x) curves, that is, dM/dx >0atg<3/2 and dM/dx < at 
q > 3/2, see Fig. 2. 




FIG. 2. Dependence of the total mass of star with first-order phase transition on the relative 
radius of new-phase core, in NNG. Abscissae are 10 6 r n /R, ordinates are 10 13 (M — 1). The upper 
panel corresponds to the case of positive a = .01. Curves (starting from upper one) are shown for 
three values of q = 1.44, 1.445 and 1.45, respectively. Note, that, in principle, at a > 0, all curves 
at x = point have negative derivative, that is there is a minor region of small-core instability for 
any q > 1. 

The reverse situation of negative a = — .01 is shown in the lower panel. Three curves (starting 
from upper one) correspond to values of q = 1.555, 1.554 and 1.553. Here, all curves at x = point 
have positive derivative, that is there is a minor region of small-core stability for any q > 1. 

The region of "abnormal" behavior of dM/dx may be very small, a value of x, at which 
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again R = 1 and M = 1 is: 



si ( 32 ° 2 ^- 1 ) V'" (21) 

. It's worth mentioning that xl is very sensitive function of q at small enough values of a. 
It seems that despite the smallness of real values of a, and therefore the smallness of the 
" abnormal" stability region, this influence of the power-law potential on stabilty of a star is 
potentially of some major interest. 



V. CONCLUSION 

The deviations from 1/r law of NTG revealing itself in changing the value power of r lead 
to a dramatic change of the stellar stabilty against the first-order phase-transition in the 
center of star, with small new-phase cores. The steeper dependence of point-mass potential 
on r leads to the absolute instability of star at smaller core dimensions, while the lesser 
(than reciprocal) dependence of the point-mass potential leads to the "absolute stability" 
of star. Though for "real" values of a this phenomenon takes place only at the very small 
dimensions of new-phase cores, nevertheless it may, in principle, serve as some means for 
proving or rejecting such form of deviations from 1/r-law. One may note, for example, 
the microcollapses of neutron star revealing itself in pulsar timing "faults" [8]. Another 
possibility may occur in the larger astronomic structures, such as IMC [9] and supermassive 
stars [3, 5]. 
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